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Abstract 


Parameter  and  state  estimation  techniques  are  discussed  for  an 
elliptic  system  arising  in  a  developmental  model  for  the  antenna  surface 
In  the  Maypole  Hoop/Column  antenna.  A  computational  algorithm  based  on 
spline  approximations  for  the  state  and  elastic  parameters  is  given  and 
numerical  results  obtained  using  this  algorithm  are  summarized. 


1.  Introduction 

This  report  summarizes  some  of  the  results  from  an  ongoing  Langley 
Research  Center  program  directed  towards  developing  parameter  estimation 
techniques  for  flexible  systems  modeled  by  partial  differential  equations 
with  an  emphasis  on  large  space  structures.  The  intent  of  the  program  is 
to  produce  general  purpose  techniques  with  a  sound  theoretical,  basis  which 
are  computationally  efficient  while  contributing  to  Langley’s  technology 
development  program  in  large  space  antennas  [i].  Of  the  many  techniques 
available  for  parameter  estimation  in  distributed  svstems  (for  example, 
see  [2])  the  spline  based  estimation  techniques  of  [3-8]  appear  well 
suited  for  large  space  structures  applications  and  are  currently  being 
developed  to  treat  this  class  of  prolioms.  Simultaneously,  an  estimation 
problem  associated  with  the  Maypole  (Hoop/Column)  antenna  [9]  is  being 
formulated  and  will  be  solved  as  part  of  the  developmental  process.  The 
next  section  of  this  paper  describes  the  Hoop/Column  antenna  and  presents 
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the  identification  problem  being  consldezed.  The  parameter  estimation 
approach  is  then  outlined  and  discussed  in  the  context  of  the  Hoop/Column 
application.  Subsequent  sections  include  mathematical  details  of  the 
antenna  application  and  numerical  results. 


II.  The  Maypole  (Hoop/Column)  Antenna 

One  of  the  planned  activities  of  NASA's  Space  Transportation 
System  is  the  placement  in  earth  orbit  of  a  variety  of  large  space  anten¬ 
nas.  Potential  large  space  missions  for  the  next  two  decades  will  require 
antennas  and  structures  ranging  from  30m  to  20km  in  size.  Applications 
Include  communications  (mobile,  trunking,  etc.),  remote  sensing  (soil 
moisture,  salinity,  etc.),  deep  space  network  (orbital  relays),  astronomy 
(x-ray  observatory,  optical  array,  radio  telescope,  very  long  baseline 
interferometry,  etc.),  energy,  and  space  platforms. 

For  the  purpose  of  technology  development,  the  NASA  Large  Scale 
Systems  Technology (LSST)  program  office  has  pinpointed  focus  missions  and 
identified  future  requirements  for  large  space  antennas  for  communications, 
earth  sensing,  and  radio  astronomy  [l].  In  this  study,  particular  empha¬ 
sis  is  placed  on  mesh  deployable  antennas  in  the  50-120  meter  diameter 
category.  Communication  satellites  of  this  size  will  require  a  pointing 
accuracy  of  0.035°  and  surface  accuracy  of  4-8mm.  One  such  antenna  is  the 
Maypole  (Hoop/Column)  antenna  shown  for  the  100m  point-design  in  Figures 
1  and  2.  This  antenna  concept  has  been  selected  by  the  LSST  office  for 
development  by  the  Harris  Corporation,  Melbourne,  Florida,  under  contract 
to  the  Langley  Research  Center  [9]. 
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FIGURE  2  Maypole  Hoop/Column  Antenna  Reflector  Surface 

The  Hoop/Coluan  antenna  consists  of  a  knitted  gold-plated  molybdenum 
wire  reflective  mesh  stretched  over  a  collapsible  hoop  that  supplies  the 
rigidity  necessary  to  maintain  a  circular  outer  shape.  The  mesh  grid  can 
be  varied  to  meet  a  given  radio  frequency  reflectivity  requirement.  The 
annular  membrane-like  reflector  surface  surrounds  a  telescoping  mast  which 
provides  anchoring  locations  for  the  mesh  center  section  (Figure  1) .  The 
mast  also  provides  anchoring  for  cables  that  connect  the  top  end  of  the 
mast  to  the  outer  hoop  and  the  bottom  end  of  the  mast  to  48  equally  spaced 
radial  graphite  cord  truss  systems  woven  through  the  mesh  surface  [9]. 
Tensions  on  the  upper  (quartz)  cables  and  outer  lower  (graphite  epoxy) 
cables  are  counter  balanced  to  provide  stiffness  to  the  hoop  structure. 

The  inner  lower  cables  produce,  through  the  truss  systems,  distributed 
surface  loading  to  control  the  shape  of  four  circular  reflective  dishes 
(Figures  1  and  2)  on  the  mesh  surface.  Flat,  conical,  parabolic,  or 
spherical  dish  surfaces  can  be  produced  using  this  cable  drawing  technique. 

After  deployment  or  after  a  long  period  of  operation,  the  reflector 
surface  may  require  adjustment.  Optical  sensors  are  to  be  located  on  the 
upper  mast  which  measure  angles  of  retroreflectlve  targets  placed  on  the 
truss  radial  cord  edges  on  the  antenna  surface.  This  information  can  then 
be  processed  using  a  ground-based  computer  to  determine  a  data  set  of 
values  of  mesh  surface  location  at  selected  target  points.  If  necessary, 
a  new  set  of  shaping  (control)  cord  tensions  can  be  fed  back  to  the 
antenna  for  adjustment. 

It  is  desirable  to  have  an  identification  procedure  which  allows  one 
to  estimate  the  antenna  mesh  shape  at  arbitrary  surface  points  and  the 
distributed  loading  from  data  set  observations.  It  can  also  be  antici¬ 
pated  that  environmental  stresses  and  the  effects  of  aging  will  alter  the 
mash  material  properties.  The  identification  procedure  must  also  allow 
one  to  address  this  issue. 

The  authors  are  currently  developing  identification  procedures  for 
use  in  distributed  parameter  models  of  the  mesh  surface.  Considering  the 
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antenna  to  be  fully  deployed  and  In  static  equilibrium,  we  are  deriving  a 
mathematical  model  which  describes  the  antenna  surface  deviation  from  a 
curved  equilibrium  configuration  (for  preliminary  findings,  see  [lO]). 
Using  a  cylindrical  coordinate  system  with  the  z  axis  along  the  mast,  it 
is  expected  that  the  resulting  model  will  entail  a  system  of  coupled 
second  order  linear  partial  differential  equations  in  two  spatial  vari¬ 
ables.  The  coefficients  of  these  equations  are  functions  of  the  material 
properties  of  the  stretched  mesh.  The  derivation  and  computer  software 
for  this  model  are  still  under  development.  In  the  meantime,  a  simpler 
developmental  (prototype)  problem  has  been  solved  which  is  descriptive  of 
the  original  problem  and  for  which  the  software  produced  will  hopefully 
be  readily  extended  for  use  in  the  more  general  case. 

For  the  developmental  problem,  the  loading  is  assumed  to  be  normal  to 
the  horizontal  plane  containing  the  hoop  rim  and  the  mesh  surface  is 
assumed  to  be  described  by  the  static  two-dimensional  stretched  membrane 
equation  [ll]  with  variable  stiffness  (elastic)  coefficients  and  appro¬ 
priate  boundary  conditions  for  the  Hoop/Column  geometry.  Mathematically, 
in  polar  coordinates,  we  have 

-7^CrECr,9)|2]--L^[E(r.8)|H]  f(r,0)  (1) 

r 

where  u(r,6)  is  the  vertical  displacement  of  the  mesh  from  the  hoop  plane, 
f(r,0)  is  the  distributed  loading  force  per  unit  area,  and  E(r, S)  >  0  is 
the  distributed  stiffness  (elastic)  coefficient  of  the  mesh  surface  (force/ 
unit  length) .  Equation  (1)  is  to  be  solved  over  the  annular  region 
ft  ■  [e,R]  x  [0,2irJ.  Appropriate  boundary  conditions  are 


u(e,8)  -  uQ 
u(R,0)  -  0 


(2) 


along  with  the  periodicity  requirement 


u(r,0)  -  u(r,2ir) 


(3) 


where  R  is  the  radius  from  the  mast  center  to  the  circular  outer  hoop,  e 
is  the  radius  from  the  mast  to  the  beginning  of  the  mesh  surface  (see 
Figure  2),  and  Uq  is  the  coordinate  at  r  *  e  of  the  mesh  surface  below  the 
outer  hoop  plane. 

We  further  assume  that  the  distributed  loading  along  with  a  data  set 

of  vertical  displacements,  u  (r.,6.),  at  selected  points  (r,,0.)  on  the 

m  j.  j  l  j 

mesh  surface  is  known.  Given  this  Information,  the  developmental  problem 
is  to  estimate  the  material  properties  of  the  mesh  as  represented  by 
E(r,0)  and  produce  estimates  of  the  surface  represented  by  u(r,0)  at 
arbitrary  (r,0)  points  within  ft.  The  procedure  applied  to  solve  this 
problem  is  discussed  in  the  next  section. 
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III.  The  Parameter  Estimation  Approximation  Scheme 

The  first  two  authors  and  their  colleagues  have  derived  techniques 
for  approximating  the  solutions  to  systems  identification  and  control 
problems  involving  delay  equation  models  and  partial  differential  equation 
models  in  one  spatial  variable  and  have  used  them  in  a  variety  of  appli¬ 
cations  [l2,  13  J.  The  Hoop /Column  application  requires  an  extension  of 
the  theory  and  numerical  algorithms  to  elliptic  distributed  systems  in 
several  spatial  variables.  The  approach,  when  specialized  to  the  param¬ 
eter  estimation  problem,  may  be  summarized  as  follows.  (1)  Select  a  dis¬ 
tributed  parameter  formulation  containing  unknown  parameters  for  a  specific 
system.  (2)  Mathematically  "project"  the  formulation  down  onto  a  finite 
dimensional  subspace  through  some  approximation  procedure  such  as  finite 
differences,  finite  elements,  etc.  (3)  Solve  the  parameter  estimation 
problem  within  the  finite  dimensional  subspace  obtaining  a  parameter  esti¬ 
mate  dependent  upon  the  order  of  the  approximation  embodied  in  the  sub¬ 
space.  (4)  Successively  increase  the  order  of  the  approximation  and,  in 
each  case,  solve  the  parameter  estimation  problem  so  as  to  construct  a 
sequence  of  parameter  estimates  ordered  with  increasing  refinement  of  the 
approximation  scheme.  (5)  Seek  a  mathematical  theory  which  provides  con¬ 
ditions  under  which  the  sequence  of  approximate  solutions  approaches  the 
distributed  solution  as  the  subspace  dimension  Increases  with  a  convergent 
underlying  sequence  of  parameter  estimates. 

In  applying  this  approach  to  the  developmental  problem  described  in 
this  report,  the  stiffness  function  is  parametrized  in  terms  of  cubic 
splines;  this  converts  the  estimation  of  E(r,0)  into  a  finite  dimensional 
parameter  estimation  problem.  After  writing  the  energy  functional  generic 
to  the  membrane  equation,  we  use  the  Galerkin  procedure  to  project  the 
distributed  formulation  onto  a  finite  dimensional  state  subspace  spanned 
by  tensor  products  of  linear  spline  functions  defined  over  ft.  The  approx¬ 
imate  displacement  thus  obtained  is  expressible  in  terms  of  the  spline 
basis  functions.  The  Galerkin  procedure  in  this  case  yields  algebraic 
equations  which  define  the  displacement  approximation  coordinates  in  terms 
of  the  unknown  parameters.  In  order  to  solve  the  approximating  parameter 
estimation  problem,  the  parameters  defining  E(r,0)  are  chosen  so  that  a 
least  squares  measure  of  the  fit  error  between  the  observed  and  predicted 
(by  the  approximate  state)  data  set  is  minimized.  Finally,  following 
steps  (4)  and  (5)  an  algorithm  is  constructed  to  determine  the  order  of 
the  linear  spline  approximation  above  which  little  or  no  further  improve¬ 
ment  is  obtained  in  the  unknown  parameters  as  one  increases  the  dimension 
of  the  subspaces. 

Following  this  procedure,  a  one-dimensional  version  of  the  develop¬ 
mental  problem  has  been  solved  [l4]  in  which  the  schemes  proposed  were 
successfully  tested.  We  have  subsequently  considered  the  two  dimensional 
case  and  details  of  these  investigations  are  presented  in  the  following 
sections. 


IV.  Finite  Dimensional  Approximations 

We  choose  a  Galerkin  procedure  [l5,  16]  with  linear  spline  basis 
functions  to  perform  the  finite  dimensional  approximation  for  the 
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developmental  problem  in  which  E(r,6)  and  u(r,9)  of  equation  (1)  are  to  be 
estimated.  The  boundary  conditions  (2)  are  first  converted  to  homogeneous 
form  by  introducing  the  new  dependent  variable 

y(r,0)  -  u(r,0)  -  (  uQ  (4) 


whereby  equation  (1)  becomes 


.  a  rE(r,9)un 


(5) 


with  boundary  conditions 
y(e,9)  -  0 

y(R,6)  -  0  (6) 

y(r,0)  -  y(r,2w)  . 


Following  the  standard  formulation  (see  [15,  16])  for  the  weak  or 
variational  form  of  (5),  the  energy  functional  £  associated  with  (S)  is 


£(*) 


2ir  R  r  .  ^ 

/  /  ^  E(r,9)7z«  7z  -  f  (r,9)z^  rdrd9 


(7) 


where  7  is  the  gradient  in  polar  coordinates  which,  in  the  form  used  here, 
is  equivalent  to 


11  J 
r  99  ' 


The  function  It  is  given  by 


(8) 


f(r,9) 


f(r,9)  +  -  T"  ( 
r  3r 


rE(r,9)uQ 
e  -  R  * 


(9) 


and  the  vertical  displacement  z(r,9)  of  the  mesh  surface  away  from  the 
hoop  equilibrium  plane  is  a  function  satisfying  the  boundary  conditions 
(6)  and  possessing  first  derivatives  on  Q  in  the  mean  square  sense  (we 

denote  this  by  z  6  hJ  _ (ft)  =  Z).  The  first  variation  6£  of  £  about  the 

u,per 

function  y(r,0)  is  given  by 


2ir  R  s  %  -s 

8E(y;v)  -  /  /  ^E-(r,9)7y  •  7v  -  f  (r,9)v|»  rdrd9 

2ir  R  r  % 

■  /  /  -^E(r,9)7y*  7v- [f  (r,9)v+  E(r,9)k  •  7v]^»  rdrdQ  (10) 
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where 


and  v  is  an  arbitrary  function  in  Z  *  HI  (0) . 

0,per 

A 

Given  a  finite  dimensional  subspace  Z  of  Zj,  the  Galerkin  procedure 
defines  the  approximation  y  as  the  solution  in  Z  of 


2ir  R  r 

I 

Os'* 


2ir  R  f 
/  /  *<  f  (r,0) 

0  e  >• 


0)7y  •  7v/*rdrd0 


v+E(r,0)k*  7v>  rdrd0  (12) 


0  e 


for  all  v  €  Z. 

For  computational  efficiency,  the  basis  functions  used  for  the  repre¬ 
sentations  of  y  in  (12)  are  taken  as  tensor  products  of  linear  B-splines 
[l5y  p.  27;  16,  p.  100].  Thus  v  and  y  are  in  the  space  spanned  by 

v“’N,  where 

v“jN(r,0)  -  o“(r)0j(0),  (i-l,...,M-l;  j-l,...,N),  (13) 

where  -  a^(r),  (i-  1, . . . ,M-1) ,  and  0^  -  0^(0),  (j  -  1, . . . ,N-1) ,  are 

standard  linear  B-splines  with  knots  uniformly  spaced  over  [e,R]  and 

[0,2ir],  respectively.  The  elements  {a . }  are  modified  to  satisfy  homo- 

N  1 

genous  boundary  conditions  while  0^  has  been  altered  to  satisfy  periodic 
boundary  conditions. 

M  N  H  N 

For  y  *  (r,0)  within  the  subspace  spanned  by  v"’  we  can  write 

U  U  M— 1  N  u  U  M  11 


M,N,  c  r  M,  .  M,NaN,fl. 

y  (r,0)  -  l  l  o  (r)w  !  6.(6)  . 

i»l  J-l  J  J 


Replacing  y(r,0)  in  (12)  by  y  *n(r,0)  from  (14)  and  successively  setting 
14  N 

v(r,0)  »  v.*  (r,0)  for  i*l,...,M-l  and  j»l,...,N  leads  to  a  set  of 
*j  M  N 

linear  algebraic  equations  for  the  w^’  coordinates  [l7].  For  large  M 

and  N,  these  equations  may  require  sparse  matrix  techniques.  We  can  avoid 
sparse  matrix  methods  by  imposing  a  separability  condition  on  E(r,0): 

E(r,0)  -  Ex(r)E2(0)  .  (15) 

M  N 

As  shown  in  [l7],  condition  (15)  reduces  the  w  *  calculation  to  the 
solution  of  the  matrix  equation  ^ 
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with 


and 


(16) 

W*. 

/  MN» 

'  (wij> 

(17) 

F- 

(/  E2(0)s”c@)sj(0)d0) 

(18) 

*- 

(/  Bl<r)Cto  °i(r)D  Cto  ap<rDrdr) 

(19) 

(/*  E2(«)[A  bJ(0)][^  eqC6Dd0) 

(20) 

/  “j(r)  0p(r)  \ 

U  Vrt  r  *) 

C 

(21) 

(  /  /R  f(r,0)a“(r)  s”(0)rdrd0 
x  n  1  1 

/2W  /  E(r,0)k6j(0 a“(r)]rdrd0)  , 


0  e 


(22) 


where.  In  (17) -(22),  i,p  -  1, . . .  ,M-1  end  j,q-l,...,N.  The  coefficient 
matrices  of  W*®*  in  (16)  have  numerically  attractive  properties: 

(i)  all  are  symetric,  (ii)  B**,  (SN,  and  D**  are  banded  (tri-diagonal),  and 

^i|  #VM 

(ill)  A  and  D  are  positive  definite. 

Research  to  construct  a  numerical  algorithm  for  solving  (16)  which 
utilizes  these  properties  is  planned.  At  present,  (16)  is  rewritten  in 
the  equivalent  form 

[(iffW®  +  W^tW1]  -  (23) 

and  solved  by  the  Bartels-Stewart  algorithm  [18]. 

In  order  to  estimate,  via  a  numerical  scheme,  the  functional  coef¬ 
ficients  E^  and  E2,  we  must  further  parametrize  these  functions  so  that 

identification  is  performed  over  a  finite- dimensional  (instead  of  an 
infinite-dimensional)  parameter  set.  To  this  end,  we  let 
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E1(r) '  Vt<r) 


(24) 


N, 


E,(0)  -  l  6.y,(e) 
j-1  3  2 


(25) 


where  v,  and  6  are  scalar  parameters  and  X-  and  \iA  are  cubic  B-spline 

K  J  R 


functions  defined  [15,  p.  6l]  over  [e,R]  and  [0,2ir^,  respectively,  whose 
orders  are  independent  of  M  and  N.  The  functions  have  been  modified  so 
that  Uj  and  its  derivatives  satisfy  periodic  boundary  conditions.  Sub¬ 
stitution  of  (24),  (25)  into  the  matrices  A^  through  E?®  greatly  simpli¬ 
fies  our  computational  efforts  since  many  components  of  these  matrices  may 
be  computed  one  time  only  in  advance  of  the  estimation  process. 

We  turn  next  to  implementation  of  the  parameter  estimation  scheme. 


i.e.,  the  numerical  determination  of  v^,  k-  1,  ...,M^,  and  5^,  j  *  1, 


N^,  that  appear  in  (24)  and  (25)  and  correspond  to  "optimal"  values  of  the 


parameters  and  E^. 


V.  The  Parameter  Estimation  Algorithm 

Appealing  to  the  ideas  found  in  previous  sections,  we  mav  now  detail 
an  algorithm  for  estimating  the  coefficients  k*l,  ...,M^,  and  6^, 

j-1,  .  ..,N_,  for  E(r,0)  that  provide  the  "best  fit"  between  approxima¬ 
tions  for  the  state  u  and  observed  data  u  obtained  from  various  sample 

IQ 

points  on  the  surface.  We  may  equivalently  consider  data  for  y  by  making 
the  transformation 

VvV  '  u»(ri'V  ‘  (  ht)uo  <26 

for  i  -  1,  ....  L  .  and  j  -  1,  . . . ,  L„. 

r  o 

We  organize  the  parameter  estimation  algorithm  into  the  following 
steps. 

1.  Select  an  order  of  approximation  for  the  cubic  spline  elements  X^, 
k- 1,  ...,M^,  and  Pj,  j-1,  ...,N^,  used  to  represent  E^  and 

Set  n  -  1. 

2.  Select  M  and  N,  the  orders  of  the  linear  spline  basis  elements  used 
to  represent  u  (and  y  ) . 

3.  As s lane  a  nominal  set  of  values  for 
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<V  V2,  ....  vM^) 


62 . * 


4.  Calculate  the  coefficient  matrices  in  (23)  and  solve  for  W^(v,6). 

M  N 

5.  Calculate,  from  (14),  y  *  v»$)  an<*  evaluate 

Lr  L9 

J^Cv.d)  -  [  [  [yM,NCri,e  ;  v,6)  -  y  (r  ,9  )]2  .  (29) 

1-1  j-1  1  j  mil 

6.  Proceed  to  step  8  If  J**’**(v,6)  is  sufficiently  small.  Otherwise, 
through  an  optimization  procedure,  determine  a  new  pair  (0,3)  which 

decreases  the  value  of  If  no  such  pair  can  be  found,  go  to 

step  8. 

7.  Set  (\>,5)  -  (v,6)  and  return  to  step  4. 

8.  Preserve  the  current  values  of  and  the  corresponding  (v,6)  pair 

as  the  n  entry  in  a  sequence  of  these  pairs,  ordered  with  increasing 
M  and  N. 

9.  Proceed  to  step  10  if  sufficient  data  has  been  obtained  to  analyze  the 
sequences.  Otherwise,  set  n  -  n+  1  and  return  to  step  2  with 
increased  M  and  N.  The  current  values  of  (v,6)  will  be  used  as 
initial  values  for  the  next  optimization  process. 

10.  From  analysis  of  the  numerical  sequences,  select  the  (M,N)  entry  which 
indicates  the  best  numerical  results.  The  corresponding  (v,S)  pair 
yields  E(r,9)  which  determines  the  material  properties  of  the  antenna 

mesh.  The  matrix  W^v  ,6),  when  used  in  conjunction  with  (14),  deter- 

M  N 

mines  an  approximation  y  '  of  the  shape  of  the  antenna  surface. 

A  convergence  theory  for  the  parameter  estimation  algorithm  may  be  found 
in  [l7].  Numerical  results  are  described  in  the  next  section. 


VI.  Numerical  Results 


Experimental  data  for  the  Hoop/Column  antenna  is  not  available  at 
this  time.  Therefore  simulated  data  is  constructed  to  evaluate  the  pre 
ceeding  algorithm. 

As  shown  in  Figure  2,  the  parent  reflector  has  four  separate  areas 
of  illumination  on  its  surface.  Each  separate  area  is  assumed  to  have 

the  same  parabolic  shape  given,  for  0  <_  0  £  and  e  <_  r  <_  R,  by 
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The  cubic  polynomial  fits  (32)  and  (34)  are  used  to  ensure  smoothness  in 

0,  in  regions  near  0  •  it,  2ir. 

It  is  expected  that  the  mesh  will  be  stiffest  near  the  outer  hoop 
(r  ■  R)  and  around  the  inner  radius  (r  ■  e) .  For  this  reason  we  choose  a 
known  value  of  E^(r)  as 

E^(r)  -  2t  -  Tsin[ir  -^j~)  ]  (e  1  r  £  R)  (37) 

A 

where  t  is  a  constant  dependent  on  the  mesh  material.  The  stiffness  in 
the  angular  direction  is  expected  to  be  uniform.  Thus  a  known  value  of 
£2(0)  is  taken  as 

E2(0)  =  T  .  (38) 

Fran  data  provided  in  [9],  a  reasonable  value  for  x  (given  in  units  /N/m  ) 
is 

t  -  3.391  ;  (39) 


similarly,  other  parameters  are  estimated  to  be  u.  ■  -7.5m,  e  ■  8.235m, 


and  R  ■  50m. 

A  10  x  24  grid  of  data  points  u^(r^,0j)  is  calculated  by  determining  u 

at  (r^»0j)  where  the  values  of  0^  correspond  to  data  taken  along  every 

other  radial  cord  truss  system;  the  10  values  of  r^  are  spaced  uniformly 

in  (e,R).  Distributed  loads  are  obtained  by  substituting  (36)-  (38)  into 
(1)  and  evaluating  for  f(r,0).  Equations  (24)  and  (.25)  are  used  with 
•  N^  ■  4  to  represent  E^(r)  and  E2(0). 

All  numerical  results  presented  were  computed  on  the  CDC  6600  at 
Southern  Methodist  University.  We  would  like  to  express  our  sincere 
appreciation  to  David  Krakosky  at  S.M.U.  who  assisted  in  the  preparation 
of  data  and  graphical  display  of  these  results.  The  optimization  scheme 

employed  to  minimize  J*®*  is  the  IMSL  version  (ZXSSQ)  of  the  Levenberg- 
Marquardt  algorithm  [l9j  ,  where  we  typically  use  default  values  of  IMSL 
parameters. 

For  each  example  reported,  values  of  N  are  specified;  in  order 
that  N  basis  elements  are  used  in  both  r  and  0  directions,  we  take 

M  ■  N+l.  Two  measures  of  performance  will  be  given  in  each  case: 


(which  provides  a  good  measure  of  state  approximation)  and  R 
using  | •  |  to  denote  the  L2  norm  on  [e,R]  *  [0,2ir], 


,M,N 


where , 


„M,N 


IE™  -  El 


|E| 


100Z 


measures  the  relative  error  between  the  "true"  parameter  E  and  the 
"optimal"  parameter  E^  associated  with  the  (M,N)C^  approximate  parameter 
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estimation  problem.  Throughout  we  choose  a  starting  guess  E®  for  the 
optimization  procedure  in  the  case  of  N-4;  for  N-8,  16,  24  and  30,  pre¬ 
vious  converged  values  of  E™  are  used  as  initial  guesses.  For  example, 
the  initial  guess  for  N*8  is  E®*  where  5f»4. 


Example  1:  He  estimate  E2(0)  only,  holding  E. (r)  fixed  at  the  true  value 
-  i  o 

E^(r).  The  starting  guess  is  EjO)  ■  1  +  .5cos6  .  Our  results  are  sum¬ 
marized  below  and  in  Figure  3. 


N 

*M,N 

a*)1* 

CP  time 

4 

5.15% 

.609 

5.41 

8 

5.72% 

.507 

5.42 

16 

5.92% 

.514 

21.58 

24 

5.95% 

.509 

55.69 

30 

5.96% 

.502 

172.99 

Example  2:  We  now  hold  £^(6)  fixed,  ££(9)  =  £2(6),  and  estimate  E^(r) 

from  the  initial  guess  of  E?  S  1.  Our  findings  are  outlined  below  and  in 
Figure  4. 


N 

jM,* 

CP  time 

4 

39.38% 

.551 

12.22 

8 

4.19% 

.418 

26.16 

16 

8.16% 

.419 

67.11 

24 

8.85% 

.411 

116.06 

30 

8.30% 

.401 

228.02 

Example  3:  We  estimate  both  E. (r)  and  E.(0).  Initial  guesses  are  given 

0  0  1  * 

by  E^(r)  ■  5  and  £2(6)  ■  1-  .  25sin0,  respectively.  In  each  case,  the 
first  coefficient  of  E2  is  held  fixed. 


N 

CP  time 

4 

42.10% 

.551 

26.88 

8 

4.21% 

.418 

36.24 

16 

8.17% 

.419 

74.23 

24 

8.87% 

.411 

190.83 

30 

8.29% 

.401 

206.02 

We  note  (by  comparing  values  of  jM,N 

for  Examples  1  and  2)  that  we 

have  greater  success  here  in  estimating  E2<0)  than  E^(r).  This  disparity 


appears  to  be  a  function  of  the  quantity  and  placement  of  sample  data. 
(See  for  example  [17]  where,  by  changing  the  data,  we  are  able  to  better 
estimate  E^  than  E2.  In  fact,  in  comparing  the  results  of  Example  2  with 

Table  6  In  [17],  it  is  easy  to  see  how  a  different  distribution  of  data 
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for  this  example  can  actually  reduce  the  values  of  shown  in  Example  2 

by  one-half.  Regardless  of  the  differences  in  estimating  E^  and  E2,  how¬ 
ever,  we  remark  that  in  all  examples  we  have  been  able  to  successfully 
estimate  the  surface  shape  of  the  model  antenna  (as  is  evidenced  by  values 

of  J*  computed  using  240  data  points) . 

We  have  also  obtained  similar  results  in  the  case  where  random  noise 
(approximately  5%  noise  level)  has  been  added  to  the  data.  These  and 
other  findings  are  sumnarized  in  Section  VI  of  [17]. 


*  —  vtttMtN  #2  <n*€> 

FIGURE  3 


FIGURE  4 
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